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ABSTRACT 

To examine the relationship between black hole mass and host galaxy velocity dispersion for low black hole 
masses, we have measured the velocity dispersions of 15 Seyfert 1 galaxies from the catalog of Greene & Ho 
(2004). These Seyferts were selected from the Sloan Digital Sky Survey to have estimated black hole masses 
below 10'' Mq. The data are consistent with a straightforward downward extrapolation of the local Mbh-o" 
relation, indicating that this correlation extends over a range of more than four orders of magnitude in black 
hole mass. The rms scatter of the sample about the extrapolated Mbh - cr relation is 0.57 dex, consistent with 
the expected scatter of single-epoch mass estimates for Seyfert 1 galaxies. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
Seyfert 



L INTRODUCTION 

The correlation between black hole mass and stellar veloc- 
ity dispersion in nearby galaxies (the Mbh - o" relation) is the 
foundation upon which our understanding of supermassive 
black hole demographics is built iF errarese & Merritt; l2000l: 
iGebhardt et ai]l2000at iTremaine etaL.2002.) . To date, how- 
ever, our knowledge of the Mbh - c relation is largely lim- 
ited to galaxies with massive bulges, velocity dispersions of 
- 70-400 km s"', and Mbh > 3 x lO*" Mq. Elucidating the 
demographics of lower-mass black holes, below 10^ Mq, can 
address key questions regarding the formation environments 
of massive black holes: for example, is there a minimum 
galaxy mass or velocity dispersion below which black holes 
are unable to form or grow ( Bromlev. Somerville. & Fabiq3 
12004)? Are black holes easily ejected from dwarf galax- 
ies by gravitational recoil following black hole mergers 
(|Merritt et al. 2004)? Extending the local black hole mass 
function to masses below 10^ Mq would also aid in predic- 
tion of event rates for black hole mergers detectable by LISA 
(e.g.,|Wvithe & Loeb 2003). 

If black holes do exist in the centers of dwarf galaxies and 
very late-type spirals, then the subset that are undergoing ac- 
cretion can potentially be detected as active galactic nuclei 
(AGNs). However, very few AGNs have been identified with 
black hole masses below IO^'Mq or in h ost galaxies with a < 
60 km S-' . These include NGC 4395 (Kraemer et alJ ll999t 
iKlippenko & Ho 2003), POX 52 (Barth et al. 2004), possibly 
NGC 4051 fcShemmer et al. 2003), and a few other narrow- 
line Seyfert 1 galaxies (Botteetal. 2004b; Grupe & Mathur 
|2004). The first systematic survey for AGN s with low-mass 
black holes was recently carried out by Greene & Hd J2004h 
(hereinafter GH), who examined the Sloan Digital Sky Survey 
(SDSS) Data Release One (DRl) archives to find 19 Seyfert 1 
galaxies having estimated black hole masses below lO*" Mq. 

The stellar velocity dispersions for the GH sample can- 
not be determined from the SDSS spectra, because the in- 
strumental dispersion of ^ 70 km s"' sets a practical lower 
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limit to the range of velocity dispersions that can be mea- 
sured, and in most cases the spectra are dominated by AGN 
emission. Here, we present new measurements of the stel- 
lar velocity dispersions for 13 galaxies from this sample and 
two additional Seyferts we identified from SDSS DR2 meet- 
ing the same selection criteria: SDSS J080907.58H-441641.4 
and SDSS J215658.30H-1 10343.1. 

2. DATA AND MEASUREMENTS 

2.1. Observations and Data Reduction 

The observations were obtained during the nights of 2003 
November 24, 2004 February 15-17, and 2004 October 9- 
10 UT at the Keck II telescope, using the ESI spectrograph 
JSheinis et al.ll2()02l) . A 0."75 slit width was used, giving an 
instrumental dispersion of ct, « 22 km s~' . The spectra cover 
the wavelength range 3800-10900 A across ten echelle or- 
ders, at a uniform scale of 11.5 km s"' pixel"'. For each ob- 
servation, the slit was oriented at the parallactic angle. Flux 
standards and red giant stars (G8III-K5III) for use as velocity 
templates were observed each night during twilight. Spectra 
were extracted with a 1" extraction width, and wavelength- 
and flux-calibrated using standard techniques. 

2.2. Stellar Velocity Dispersions 

We measured the stellar kinematics by fitting broadened 
stellar spec tra to the galaxy spe ctra, following techniques de- 
scribed by Barth. Ho. & Sargent (,2002). A Gaussian veloc- 
ity broadening function was assumed, and the fits included 
both an additive featureless continuum component and a mul- 
tiplicative low-order polynomial that can compensate for red- 
dening in the galaxy spectrum. Fits were performed in two 
spectral regions: ~ 5050 - 5400 A, which includes Mgb, 
Fe5270, and numerous weaker features; and ^ 8450-8700 
A, covering the Ca II near-infrared triplet. The exact fitting re- 
gion was tailored for each galaxy to fit within a single echelle 
order and to exclude emission lines such as [N I] A5200 and 
[Fe VII] A5158. The measurement uncertainty was taken to 
be the sum in quadrature of the fitting uncertainty from the 
best-fitting template and the standard deviation of results from 
all templates. A minimum of six template stars were used for 
each measurement. We note that the problem of template mis- 
match in the Mg/Fe abunda nce ratio, which is often n oticeable 
in high-dispersion galaxies dBarth. Ho. & Sargentl2002j) . does 
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TABLE 1 

Observations and Measurements 



— 

Object 




— 

Exposure 

(s) 


S/N 


— 7'- — 7. ; — 

Line Strength 


(km s"' ) 


,^ — 

cr(Ca 11) 

(km s"' ) 


rVvriM(rlQf) 
(km s"' ) 


(5100 A) 


log(MBH/M0) 


GHOl 


077 


2982 


17 


n 7A-|-n OR 


36 ± 6 


34 ± 8 


950 


43 1 


6.28 


GH02 


030 


6000 


21 




50 ± 4 


59 ± 8 


690 


41 9 


5.16 


GH03 


102 


3000 


10 


n 45 -1- n 07 

Lf.t J yj.\f 1 


55 ± 6 




820 


42 9 


6.01 


GH04 


043 


2700 


32 


n 77 _i_ fi 05 


(^0 -1- 4 


57-1-7 




42 8 


5 87 


GH05 


074 


3600 


26 








720 


43 1 


6 04 


GH06 


0.100 


1800 


10 


0.36 ±0.06 


48 ±7 




1240 


42.9 


6.37 


GH07 


0.094 


3600 


17 


0.17 ±0.04 


39±6 




1140 


43.0 


6.37 


GHIO 


0.081 


1800 


29 


0.17 ±0.04 


49±7 


54±8 


720 


43.3 


6.18 


GHll 


0.081 


1800 


19 


0.37 ±0.08 


64±6 


62±7 


2110 


43.2 


7.04 


GH13 


0.126 


1800 


12 








980 


43.3 


6.45 


GH14 


0.028 


1800 


14 


0.34 ± 0.05 


48 ±7 


59±7 


770 


41.9 


5.26 


GH16 


0.069 


1100 


19 


0.32 ±0.07 


69±6 


83 ±8 


1170 


42.5 


6.04 


GH17 


0.099 


3500 


20 


0.27 ±0.06 


58±7 




930 


43.1 


6.26 


GH18 


0.183 


3600 


12 


0.49 ±0.17 


81±9 




1770 


43.1 


6.81 


GH19 


0.036 


3600 


20 


0.69 ±0.19 


58±4 


53±4 


1890 


41.7 


5.89 


0809+4416 


0.054 


3600 


22 


0.47 ±0.11 


65±4 




1100 


42.8 


6.19 


2156+1103 


0.108 


3600 


34 


0.23 ±0.06 


81±8 




780 


43.5 


6.39 



Note. — Object names are those given in the Greene & Ho ( 2004) catalog. S/N is the mean signal-to-noise ratio per pixel in the ESI spectrum in the 
spectral region around Mgfo used to measure cr. "Line Strength" gives the strength of stellar absorption features in the AGN spectrum relative to those in 
the template stars, over the blue fit region. FWHM(Ha) is the width of the broad component of Ha, as determined by multi-Gaussian model fits. The AGN 
continuum luminosity log[AL;^(5100 A)] is from lGreene & Hcl 120041) : \Lx is in units of erg s"'. The luminosities are corrected for Galactic extinction 
and are calculated for a cosmology with Hq = 72 km s"' Mpc"' , f2m = 0.3, and Qa = 0.7. 



not significantly affect this sample. 

We were able to determine velocity dispersions for 13 of 
the 15 objects we observed from the GH sample, as well as 
the two new SDSS objects, using the Mgb spectral region. 
The other two galaxies were dominated by AGN emission 
and no useful results could be recovered from the template 
fitting. For 7 galaxies we were able to measure a from the 
Ca II triplet region as well; the results are generally consistent 
with the M.gb results but have larger uncertainties because of 
the smaller spectral region used for the fitting. Thus, we use 
the Mgb measurements as the best estimates of the stellar ve- 
locity dispersions. The fitting results are listed in Tableland 
illustrated in Figure^ (The figure displays, for clarity, only a 
portion of the fitting region.) 

We are unable to determine whether the measured velocity 
dispersions are dominated by light from bulges, disks, or nu- 
clear star clusters, since in most cases the galaxies are only 
marginally resolved in the SDSS images. High-resolution 
imaging with the Hubble Space Telescope (HST) would clar- 
ify the host galaxy morphologies and would also make it pos- 
sible to examine the behavior of the Mbh - i^buige relation for 
this sample. 

2.3. Black Hole Masses 

Black hole masses were derived by GH from the SDSS 
spectra using the empirical relations from Kaspi et al. ( 2000). 
These relations combine the broad-line region (BLR) radius 
(derived from the nonstellar 5100 A continuum luminosity) 
and broad-line width to obtain a virial estimate of the cen- 
tral mass. Single-epoch virial estimates of Mbh based on this 
technique are thought to have a random scatter of roughly a 
factor of 3, based on the consistency of the masses with the 
Mbh-o- relation and other tests (e.g., Gebhardt et al.' 2000b"; 
Eerrarese et al. 2001; Vestergaard 2002; Nelson et al. 2004; 
[Onken et al. 2004). While FWHM(H/5) is most commonly 
used as a measure of the velocity distribution of the BLR, the 
H/3 emission lines in the SDSS spectra have low S/N and GH 



used FWHM(Hq!) as a surrogate. 

The high S/N and spectral resolution of the ESI data make 
it possible to perform decompositions of the Hq;H-[N II] lines 
that are more accurate than the previous measurements from 
SDSS spectra. Following the methods used by GH, we fit 
Hq;H-[N II] using a multi-Gaussian model. All three narrow 
components (Ha and [N II] AA6548,6583) were constrained 
to have the same velocity width, and the wavelength separa- 
tion and intensity ratio of the [N II] lines were fixed at their 
laboratory values. Two Gaussians were required to fit the 
broad component of Ha adequately for most galaxies. In 
most cases our results for FWHM(Hq;) are similar to those 
of GH, but in a few cases such as GHll the ESI spectrum 
yields a substantially larger width, which we attribute to the 
difficulty of deblending the line profiles in SDSS spectra hav- 
ing weak broad components and only marginally resolved nar- 
row lines. The difference in linewidths could also be due in 
part to intrinsic variability. The statistical uncertainties on 
FWHM(Hq;) from the profile fits are small, typically only 
^ 5 -20 km s~', while the true uncertainties on the linewidths 
are dominated by systematics including starlight subtraction 
errors, and we estimate these uncertainties to be of order 10% 
in FWHM(Hq!). We also measured the H/3 linewidths using 
a double-Gaussian model to represent the broad component. 
For three objects the broad H(3 line did not have sufficient S/N 
to obtain a meaningful fit, and in one object (GHl 1) broad H/3 
was not detected at all. We do not find any systematic trend 
for H/3 to be significantly broader than Ha, and we use the 
Ha lines to determine Mbh because their higher S/N enabled 
us to measure the linewidths in a consistent manner for the 
entire sample. 

Since the Keck spectra were observed and extracted 
through a small aperture and the slit losses are unknown, 
we continue to use the GH measurements of the AGN con- 
tinuum luminosity obtained from the SDSS spectra. These 
may slightly overestimate the true AGN luminosity because of 
residual starlight in the SDSS spectra. HST V-band imaging 
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Fig. 1 . — ESI observations of the spectral region around Mgb and Fe5270. 
The spectra are flux-calibrated in fx units and have been normalized to a 
median flux level of unity in order to show the varying amounts of nonstellar 
continuum dilution in the sample. The red overplotted spectra are the best- 
fitting broadened template stars. The spectrum of the K2111 star HD 332344 
is shown at the bottom right; note the expanded vertical scale of this panel. 



of this sample would be the best way to obtain more accurate 
luminosities. 

Single-epoch virial mass estimates for black holes in AGNs 
are subject to potentially significant systematics, due to un- 
certainties in both the slope 7 of the BLR radius-luminosity 
relation (rsLR oc L^) and the normalization factor / in the 
virial relation Mbh = /'"blr^fwhm/^- The radius-luminosity 
relation is not well determined for XL\ < 10'^^-^ erg s"', and 
this presents particular problems for our sample since it is 
necessary to extrapolate this relation to luminosities as low 
as ALa < 10"^' ^ erg s"'. Following Kaspi et al. (2000) and 
IVesterg aard (2002) we use 7 = 0.7; setting 7 = 0.5 instead 
(e.g.. [Shields et al. IQ Q'},) would raise the masses by 0.3 dex 
on average. To obt ain virial mas s es, we use the recent cal- 
ibration of / from lOnken et al.l (120041) and iPeterson et al.l 
(|2004), which was derived by scahng the black hole masses 
of reverberation-mapped Seyferts to fit the Mbh ~ cr relation. 
This calibration increases the black hole masses by 0.27 dex 
relative to masses derived for the assum ption of isotropic or- 
bits of BLR clouds ('e.g.. iKaspi et al J 12000) . and shifts the 
sample to systematically higher masses than those originally 
given by GH. The revised masses are listed in Table 1 . 



3. DISCUSSION 

Figure|2ldisplays the Mbh-c relation for AGNs with low 
black hole masses as well as for nearby galaxies with di- 
rect dynamical measurements of Mbh- For the SDSS sample, 
NGC 4395, and POX 52, the vertical error bars represent a 
factor of 3 uncertainty in Mbh, corresponding to the expected 
level of random scatter for single-epoch measurements. Over- 
all, the SDSS AGNs lie close to the Mbh - cr relation extrapo- 
lated to low masses. Compared to the Tremaine et al. (2003) 
Mbh-ct relation with power-law slope 4.02, the SDSS galax- 
ies have a mean offset of 0.23 dex in Mbh and an rms scatter 



of 0.56 dex. The mean offset and rms scatter are 0.64 dex 
and 0.85 dex, respectively, relative to the steeper Mbh - c re- 
lation of ' Merritt & Ferraresd ( 120011) with slope 4.72. (The 
mean offset from the Tremaine et al. Mbh - cr relation would 
be -0.04 dex if we used the virial relations from Kaspi et al. 
2000 rather than the Onken et al. calibration.) The scatter 
of 0.56 dex relative to the shallower Mbh - cr relation is only 
slightly larger than the scatter of 0.4-0.5 dex for AGNs with 
reverberation-based black hole masses (Onken etal. 200l 
iNelson et al 2004). This overall consistency suggests that 
both the BLR radius-luminosity relation and the Mbh -cr rela- 
tion can reasonably be extrapolated to black holes in the range 
- 10^ - lO*" Mq, well below the mass scales in which these 
correlations have been calibrated. 

Given the systematic uncertainties involved in deriving 
single-epoch virial masses for these objects and the substan- 
tial scatter, we do not attempt to fit the slope of the Mbh -cr re- 
lation based on the new data. Nevertheless, our results may be 
useful as a constraint on models making specific predictions 
for the low-mass end of the Mbh - cr relation. For example, 
the hierarchical black hole growth scenario of Granato et ay 
(|2p04) predicts that the Mbh - cr slope should substantially 
steepen below cr = 150 km s"' as supernova feedback becomes 
more efficient at slowing the AGN fueling rate. The low-mass 
AGN sample does not show such behavior. 

We cannot rule out the possibility that the GH sample is 
only the upper envelope of the population of black holes in 
low-mass galaxies. The selection technique used by GH is 
only possible for broad-lined AGNs, and the probability of 
detecting a galaxy as an AGN in the SDSS archive depends 
on the AGN luminosity, which in turn depends on Mbh- This 
potentially biases the sample towards preferential detection of 
AGNs with the highest-mass black holes for a given cr. There 
are still few direct constraints on the masses of black holes 
in inactive dwarf ellipticals or late-type spirals, and the case 
of M33 demonstrates that some late-type galaxies may not 
contain a black hole at all (Gebhardt et al. 2001). Dynam- 
ical searches for low-mass black holes should be continued 
whenever feasible in order to obtain an unbiased sample of 
detections or upper limits in this mass range. 

Based on their broad-line widths of < 2000 km s"', the 
GH objects can be classified as narrow-line Seyfert 1 galax- 
ies (NLSls). Our observations provide new data to test re- 
cent claims regarding the MBH-cr relation of NLSls. Us- 
ing [O III] emission-line velocity dispersions as a surrogate 
for stellar velocity dispersions, Grupe & Mathur (2004) found 
that ROSAT-selected NLSls lie below the Mbh-ct relation, 
typically by an order of magnitude, and concluded that NLS 1 s 
are "young" AGNs with undermassive black holes. We find, 
on the contrary, that the low-mass NLS 1 galaxies in the GH 
sample do not deviate significantly from the local Mbh - cr re- 
lation. The use of [O III] emission-line velocity dispersions 
by Grupe & Mathur (2004) in place of stellar velocity disper- 
sions may be the cause of this discrepancy. As recently shown 
by Botte et al. (2004a), NLSls apparently deviate from the 
Nelson & Whittle ( 1996) empirical trend for Seyferts to have 
similar stellar and emission-line velocity dispersions. We will 
examine the relation between the stellar and emission-line ve- 
locity dispersions for the GH sample in a future paper, in the 
context of a larger sample of Seyfert galaxies covering a wide 
range of black hole masses. 

4. CONCLUSIONS 
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Fig. 2. — The Mbh ~o" relation for massive blaclc holes in galaxy nuclei. 
P oints in blac k are black holes with dynamical mass measurements compiled 
bv lTremaine et alJi2002) . Open circles represent stellar-dynamical measure- 
ments, filled triangles are gas-dynamical detections, and crosses are from 
H2O maser observations. Red filled circles are the SDSS Seyfert 1 galaxies 
describe d in this paper. F illed squares are measurements or upper limits fo r 
IC 342 ( Boker. van der Marel. & Vacca 1999), POX 52 iBarth et al. 200?), 
NGC 4395 (FiliDDenko & Ho 2003), and M33 (Gebhardt et al. 2001). For 
the SDSS sample, POX 52, and NGC 4395, the vertical error bars correspond 
to a factor of 3 uncertainty in mass. The mass for POX 52 has been scaled 
upward by 0.27 dex from the value given by Barth et al. ( 2004) to reflect the 
new virial mass calibration from Onken et al. i2004). The solid and da shed 
lines represent the Mrh - cr relations derived bv Tremaine et all 12003) and 
IMerritt & Ferrares3 1200 II) . respectively. 



We have measured the stellar velocity dispersions of 15 
Seyfert 1 galaxies selected from the SDSS to have low-mass 
black holes based on their broad-line widths and nonstellar 
continuum luminosities. The velocity dispersions for this 
sample are small, ranging from 36 to 8 1 km s"' , and the re- 
sults are consistent with a simple extrapolation of the M^^-a 
relation to low masses. This initial view of the demographics 
of low-mass black holes in AGNs implies a remarkable conti- 
nuity in the growth and evolution of black holes and their host 
galaxies for Mbh ranging from ^ 10^ to over 10^ Mq. Fur- 
ther refinement of the black hole mass estimates is possible 
if reverberation mapping can be successfully carried out for 
AGNs having Mbh < 10^ Mq\ this may be the best prospect 
for a definitive measurement of the Mbh - cr slope for low- 
mass black holes. 
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